Interplay of Aharonov-Bohm and Berry phases in gate-defined graphene quantum dots 
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We study the influence of a magnetic flux tube on the possibility to electrostatically confine elec- 
trons in a graphene quantum dot. Without magnetic flux tube, the graphene pseudospin is responsi- 
ble for a quantization of the total angular momentum to half-integer values. On the other hand, with 
a flux tube containing half a flux quantum, the Aharonov-Bohm phase and Berry phase precisely 
cancel, and we find a state at zero angular momentum that cannot be confined electrostaticaUy. In 
this case, true bound states only exist in regular geometries for which states without zero-angular- 
momentum component exist, while non-integrable geometries lack confinement. We support these 
arguments with a calculation of the two-terminal conductance of a gate-defined graphene quantum 
dot, which shows resonances for a disc-shaped geometry and for a stadium-shaped geometry without 
fiux tube, but no resonances for a stadium-shaped quantum dot with a 7r-flux tube. 

PACS numbers: 73.63.Kv, 73.22. Pr 



I. INTRODUCTION 



In recent years, graphene has emerged as a promising 
material for future nanoelectronical devices.'^"^ The pos- 
sibility to confine electrons is of particular relevance in 
this context. Experimental activity concentrates on con- 
finement in quantum dots realized with etched graphene 
~ Hi 



structures'"^ or graphene nanoflakes.'^ Electrostatic con- 
finement with the help of metal gates, which is standard 
in semiconductor heterostructures, is problematic due to 
the absence of a band gap in monolayer graphene. In par- 
ticular, an electron that approaches a region of graphene 
with zero carrier density — the closest approximation to 
an "electrostatic barrier" in graphene — will penetrate 
this region with unit probability if at normal incidence. 
This phenomenon is known as "Klein tunneling" .EHUl 
Theoretical proposals suggest to use magnetic instead of 
electric fields to shape quantum dots^^ or induce a gap 
in the spectrum.li^ 

The statement that one cannot confine electrons in 
graphene using gate potentials can be circumvented in 
certain special cases.*^ The reason is that Klein tunnel- 
ing is effective only at perpendicular incidence, while the 
reflection probability sharply increases away from nor- 
mal i nciden ce. Certain integrable geometries, such as 
a disc) ^^ * ^^ -^ allow states that exclude perpendicular inci- 
dence, so that electrons can be effectively confined in a 
disc-shaped region of graphene with finite carrier den- 
sity, surrounded by a carrier-free {i.e., undoped, intrin- 
sic) graphene sheet. On the other hand, for geometries 
with a chaotic classical dynamics, no such exclusion of 
perpendicular incidence is possible, and one may expect 
that no bound states exist in this case. In Ref. T^, as well 
as in later studies, "'^'^ ^^ ^^ a circular and a stadium-shaped 
quantum dot, as prototypes of integrable and chaotic ge- 
ometries, were embedded in a carrier-free graphene re- 
gion and coupled to source and drain contacts, as shown 
schematically in Fig. [T] Bound states are then revealed 
as sharp resonances in the two-terminal conductance. 
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FIG. 1: (Color online) The geometry under consideration: 
A quantum dot (here with circular shape), consisting of a 
region of graphene with a nonzero spatially uniform carrier 
density surrounded by a carrier-free graphene layer, which 
is coupled to leads in a two-terminal geometry. The total 
system has rectangular shape of dimension L x W; the size 
of the quantum dot is denoted R. In this article we consider 
the effect of a magnetic flux tube through the quantum dot 
(indicated in red). If the flux tube carries half a flux quantum, 
the Aharonov-Bohm phase precisely cancels the Berry phase 
that is accumulated in a cyclic orbit inside the quantum dot. 



Interestingly, the conductance of a carrier-free 
graphene sheet with a stadium-shaped and disc-shaped 
quantum dot showed resonant features t hat we re quanti- 
tatively different, but qualitatively similar .tHEl The quan- 
titative difference concerns the scaling of the resonance 
widths with the coupling to the leads, which is deter- 
mined by the ratio R/L of the quantum dot size R and 
the distance L between the source and drain contacts, see 
Fig. [1] Whereas for the stadium-shaped quantum dot 
the width was proportional to R/L for all resonances, 
the disc-shaped quantum dot also featured much nar- 
rower resonances, with a width that scaled proportional 
to (R/L)"' with n > 3. The qualitative similarity was 
that in the limit R/L — > both systems showed con- 
ductance resonances at all. This contradicts the naive 



classical expectation that there should be no resonances 
for a chaotic geometry, because in a chaotic geometry 
each electron eventually hits the dot boundary at per- 
pendicular incidence, and an electron that hits the dot 
boundary at perpendicular incidence exists the quantum 
dot with unit probability. No resonant structures should 
exist if the escape probability is unity after a finite time. 

This deviation from the naive classical expectation 
can be attributed to the Berry phase in graphene. In 
graphene, electrons are assigned a pseudospin that cor- 
responds to the sublattice degree of freedom. The pseu- 
dospin is locked to momentum. Upon completion of a 
full rotation the electron collects a Berry phase of tt. 
This Berry phase has the important consequence that 
the lowest possible angular momentum is |. Perpendic- 
ular incidence on the boundary of the quantum dot cor- 
responds to zero angular momentum, so that no states 
with perpendicular incidence on the surface exist. This 
then explains why a stadium-shaped quantum dot still 
shows conductance resonances, in spite of the naive clas- 
sical expectation that geometries with a chaotic classical 
dynamics can not be used to confine electrons. 

In order to support these arguments, in this article 
we study gate-defined graphene quantum dots in which 
the Berry phase is compensated by the Aharonov-Bohm 
phase from a magnetic flux tube through the quantum 
dot. (The interplay of the two phase shifts is also used 
experimentally to identify Berry phase effects, see, e.3., 
Refs. 18 19.) If the magnetic flux tube carries half a flux 
quantum ("tt flux"), the Aharonov-Bohm phase and the 
Berry phase collected along a closed trajectory around 
the flux tube precisely cancel. We find that with a tt flux 
the system can reach a state with zero kinematic angular 
momentum, that cannot be confined by means of gate 
potentials. 

The TT-fiux tube has qualitatively different conse- 
quences for disc-shaped and stadium-shaped geometries. 
For the disc-shaped geometry, states have a well-defined 
kinematic angular momentum. While the states with 
zero kinematic angular momentum are no longer con- 
fined, states with nonzero kinematic angular momentum 
remain confined to the quantum dot. Hence, for the disc- 
shaped quantum dot the inclusion of the 7r-flux tube elim- 
inates some of the resonances, but not all. On the other 
hand, for the stadium-shaped geometry, all states have 
a component in the zero-angular-momentum channel, so 
that inclusion of the 7r-flux tube leads to the suppression 
of all resonances.'221 

The remaining part of the paper is organized as follows: 
In Section HI] we calculate the bound states of a disc- 
shaped quantum dot in the presence of a magnetic flux 
tube. We flnd, that the asymptotic behavior of the zero- 
angular-momentum state is the same as for a free circular 
wave. Hence, no bound state can exist in this channel. In 



Section HI we present a numerical calculation of the two- 
terminal conductance setup of Fig. [IJ for a circular and 
a stadium-shaped quantum dot. Upon inclusion of the 
TT-flux tube, we find that sharp resonances persist for the 



circular dot, while the conductance becomes featureless 
for the stadium dot in the limit R/L — >■ 0. We conclude 
in Section HVl 



II. DISC-SHAPED QUANTUM DOT 

The electrostatically-defined graphene quantum dot is 
described by the Hamiltonian 



Hq = wf(p + eA) • (t + U(r), 



(1) 



where v-p is the Fermi velocity wp and cr = {ax,<^y) are 
the Pauli matrices. The gate potential V{r) is nonzero 
and constant inside the quantum dot, and zero elsewhere. 



U(r) 



— fivpUg, r < R 
0, r>R, 



(2) 



where R is the radius of the disc-shaped dot. The con- 
stant Vb has the dimension of inverse length. We choose 
Vq > 0, so that the quantum dot is electron doped. The 
potential V{r) is smooth on the scale of the lattice con- 
stant, justifying our description in terms of a single Dirac 
point. The choice of a spatially uniform potential inside 
dot makes a closed-form solution of the wavefunctions 
possible and allows for a straightforward comparison to 
the classical dynamics in the quantum d ot, bu t it is not 
essential for the existence of bound states.l^ilSll^jjg struc- 
ture of quasibound states in the inverted setup (zero po- 
tential inside, nonzero outside) was considered in Refs. 

The vector potential corresponding to the magnetic 
fiux line is 



A(r) 



h $ . 
< 

e 27rr 



(3) 



where eg is the unit vector for the azimuthal angle, and <& 
is the magnetic flux measured in units of the flux quan- 
tum h/e. In polar coordinates, the kinetic part of the 
Hamiltonian then reads 



^f(p 



= —ihv\i 











where we defined the operators 



D+ = 



,±i9 



dr ± -de T 



$ 



(4) 



(5) 



With our choice of the vector potential, the Hamil- 
tonian is invariant under rotation, hence we can look for 
eigenstates of the total angular momentum j^ = lz + ^<^z- 
They have the form 



V'm(r) 



AmO 



e '2¥?m,+(r) 



(6) 



where m = ±1/2, ±3/2, .... Inside the dot, for r < R, 
the radial wave functions (pm,± are determined by the 
coupled equations 

{dr - (to- i)i - i$) tp,n,+ ir) = iVoipm,-ir), 

{dr + ("1 + i)i + i$) (^,„,-(r) = iVo^,n,+ (r). (7) 

Outside the dot the equations decouple, and the radial 
wave functions show a power law behavior 



with coefficients a±. 



-m-l/2-* 



(8) 



A. Without flux tube 

We first review the solutions without flux tube, for 
$ — O}^ With the requirement that the wavefunction is 
regular for r ^- 0, we find for the solution inside the dot 



fm,+ (r) = J\m-l/2\{VQr), 
fni-{r) = isgn(TO)J|„+i/2| (Vbr), 



(9) 



where Jn{x) is the Besscl function. Outside the dot, the 
wave function must not diverge, which gives the con- 
straints a+ = (to > 0) and a^ — (to < 0). From 
continuity of the wavefunction at r = R, we find the 
resonance condition 



J\m\-l/2iyoR) = 0- 



(10) 



The wavefunction outside the dot is decaying as ex 

r-(l™l + l/2). 



In Section 111 we connect the quantum dots and the 



surrounding undoped graphene layer to source and drain 
contacts. The distance between the contacts is denoted 
L and the quantum dot is placed halfway between the 
contacts, see Fig. [ij In the limit L ^ R, the bound 
states are then revealed as resonances in the two-terminal 
conductance as a function of the gate potential Vq. These 
resonances have a finite width F, which can be estimated 
asPri? - \'ilj{L)\^L/\i:{R)\^R. We conclude, that the 
width of the resonances without flux tube scales as 



Fi?oc 



2|m| 



(11) 



For I TO, I = 1/2 the wavefunction decays proportional to 
1/r, which is marginally non-normalizable. Despite the 
absence of a bound state in the strict sense, the con- 
ductance nevertheless shows a resonance, with a width 
rE(x(i?/L).lSIIlII3 



B. With flux tube 

We now consider a disc-shaped quantum dot with a 
flux tube carrying half a flux quantum ($ = 1/2) — a 



"tt flux" — at its center. The results take a form similar 
to those without flux tube if we consider the kinematical 
orbital angular momentum. 



lz,kin = [i" X (p + eA)] 



(12) 



instead of the canonical angular momentum. With the 



inclusion of a 7r-flux, we then find I 



^,kin 



lz + 



f. The 



wavefunctions from Eq. (ro| are then eigenstates oi jz, kin 
with eigenvalue /ifi, where fi — m + 1/2, i.e. the kine- 
matical angular momentum takes on integer values. For 
/i ^ the calculation for the bound states proceeds in the 
same way as without flux, and we find that the resonance 
condition is given by 



JM^i/2iVoR) - 0. 



(13) 



Outside the dot, the wavefunction decays proportional 
to r~'l'^l"'"^/^). We conclude that, if the dot and the sur- 
rounding undoped graphene layer are contacted to source 
and drain reservoirs, the width F of the resonances in the 
two-terminal conductance scales as 



FiT!cx 



2|mI 



(14) 



The state with zero kinematical angular momentum 
(^ ~ 0) however is special: First of all, inside the dot, 
the wavefunction is of the form 



V'(r) = bi 



^Y,/2{Vor) 



+ &2 



-iJi/2(Vbr) 



(15) 

Recalling that the half-inte ger Be ssel functions take the 
simple form Ji/2{x) = •\/2/7rxsinx, and Yi/2{x) = 
we see that ip{r) diverges as 1/y/Vor at 



— -\/2/7rxcosa;, 
the origin, and that there is no non-trivial choice of co- 
efficients &i and 62 which removes this divergence. The 
root of this singular behavior lies in the vector potential, 
which is singular upon approaching the origin. The prob- 
lem can be cured by regularizing the vector potential. 
One possibility is to let the flux $ have an r-dependence, 
such that $ = for r < p and $ = 1/2 for r > p, i.e., 
the flux is not located at the origin, but on a circle of 
radius p. Obviously, the problem is now well-deflned at 
the origin, and we can take the solution from the case 
without flux tube, 



^(r) 



-'Vi(Vbr) 
-iJoiVor) 



(16) 



where c is a complex constant. We then match the wave- 
functions from Eq 
taking p 



(16) and Eq. (15) at r = p 



Upon 
as a condition for Eq. ([l5|. 



0, we get 62 

The boundary condition at the origin ensures, that there 
is precisely one solution for zero angular momentum. 

The p = state is also special outside the dot, where 
the wavefunction is proportional to -1= in both compo- 
nents. Thus it has the same decay as a free circular wave 
in two dimensions and, hence, it does not allow for the 



formation of a bound state. This conclusion is indepen- 
dent of the choice of the regularization of the wavefunc- 
tion near r = 0. 

Summarizing: Without flux tube, the bound states are 
labeled by the angular momentum quantum number m, 
which takes half-integer values. For |r7i| ~ 1/2 one has 
a "quasi-bound state" , because the corresponding wave- 
function is marginally non-normalizable. With a tt flux 
tube, the bound states are labeled by the kinematic angu- 
lar momentum quantum number fj,, which takes integer 
values. There is no bound state for /x = 0. 



III. TWO-TERMINAL CONDUCTANCE 



Following Refs. ll4ll6ll7l we now attach metallic source 
and drain contacts to the undoped graphene layer that 
surrounds the quantum dot. Schematically, this setup 
is shown in Fig. [l] We then calculate the two-terminal 
conductance, where bound states of the dot show up as 
resonant features as a function of the gate voltage Vq. 

The contacts are included by the addition of an addi- 
tional potential C/icads wittPS 



^ leads 



OO 



if -L/2 <x< L/2, 

if X < -L/2 or a; > L/2. 



(17) 



We apply periodic boundary conditions in the y direction, 
with period W. For the vector potential A we take a 
different gauge than in Sec. [llj 



A(r) 



/i$ 



Six) 



if < y < W/2, 

1 if -W/2 <y <0, 



(18) 



where e^; is the unit vector in the x direction. With this 
choice of the vector potential there are two flux tubes: 
one, at y = 0, located in the quantum dot, and one, at 
y = W/2, located outside the quantum dot. The second 
flux tube is necessary to implement the periodic bound- 
ary conditions. It does not affect the conductance res- 
onances in the limit that the sample width W is much 
larger than the distance L between source and drain con- 
tacts. 

The numerical calculation of the two-terminal conduc- 
tance follows the method of Ref. [2H1 Details specific to 
the presence of the flux tube are discussed in the ap- 
pendix. We now compare results for quantum dots with 
and without flux tube. We give results for a disc-shaped 
quantum dot, as a prototype of a quantum dot with in- 
tegrable dynamics, and a stadium-shaped quantum dot, 
the prototype of a dot with chaotic dynamics. 



A. Disc-shaped dot 

The two-terminal conductance for the case of a disc- 
shaped quantum dot without and with flux tube is shown 
in Fig. [2j The figure shows pronounced resonances as a 
function of the gate voltage Vq, with positions that agree 




FIG. 2: (Color online) Two-terminal conductance of a 
graphene sheet containing a disc-shaped quantum dot with- 
out (top) and with (bottom) vr-flux tube. Model parameters 
are R/L = 0.2 and W/L = 6. Without flux tube, resonances 
have definite angular momentum, with quantum number \m\ 
indicated at each resonance [data taken from Ref. 14^. With- 
out fiux tube, resonances are labeled by the kinematic angular 
momentum quantum number |/i|. No resonance is found for 
M = 0. 



with the ones calculated Sec. [TTj Without fiux tube, the 
resonances are labeled by the quantum number \m\ = 
1/2, 3/2, 5/2, .... Their width scales oc (i?/L)2|™l as 
the coupling to the leads is decreased (data not shown), 
consistent with Eq. (11) and Refs. |14| il6j j l7 ( With flux 



tube, the resonances are labeled by the kinematic angular 
momentum quantum number |/z| = 1, 2, 3, ... . There are 
no resonances for /i = 0. Upon decreasing the coupling 
to the leads, the resonances become narrower but retain 
their height, see Fig. |3J and the scaling of the resonance 
width with the ratio R/L is consistent with Eq. ( 14 ) (data 
not shown) . 



B. Stadium-shaped dot 

As a prototypical example of a chaotic quantum dot, 
we consider a stadium-shaped quantum dot. Here the 
potential T^(r) = —hvpVo for positions r inside the sta- 
dium and V{r) = otherwise. Without magnetic flux, 
the two-terminal conductance shows resonances, which, 
in the limit of small R/L, all behave as the |r7i| — 1/2- 
type resonances of the disc-shaped dot, i.e., their height 
remains finite, whereas the resonance width scales pro- 




FIG. 3: (Color online) First two resonances for a disc-shaped 
quantum dot with yr-flux tube, for different coupling strengths 
to the leads. Calculations are performed for W/L = 8 and 
various R/L, as indicated in the figure. The second resonance 
is shown enlarged in the inset. 



portional to R/LW^ The numerical data shown in the 
top panels of Figs. |4] and [5] clearly reveal these reso- 
nances, although the asymptotic scaling of the resonance 
width and resonance height with R/L is somewhat ob- 
scured by transient contributions for moderate R/L that 
originate from higher-angular-momentum contributions 
to the resonances.l^ 

The conductance trace for a stadium-shaped quantum 
dot with a flux tube carrying half a flux quantum is shown 
in the bottom panels of Figs. [4] and [5j fn order to break 
inversion symmetry, the stadium is placed asymmetri- 
cally with respect to the flux tube, see the inset of Fig.|4] 
The differences with the case of the disc-shaped quantum 
dot and with the case without a flux tube are significant. 
We find that the conductance depends on the gate volt- 
age Vo for finite R/L, but the widths of the "resonances" 
is independent of the coupling to the leads, which is set 
by the ratio R/L, whereas the height decreases upon de- 
creasing R/L. This agrees with the expectation that, 
since all states in the stadium have a /x = compo- 
nent, a stadium dot should not support any (quasi)bound 
states. While for intermediate values of R/L contribu- 
tions from higher angular momentum channels still give 
rise to broad "quasi-resonances" , in the limit R/L — >■ 0, 
only the ^ = channel is relevant, and the conductance 
becomes featureless as a function of Vq. 

We remark that, if the flux tube would be placed ex- 
actly in the middle of the stadium, inversion symme- 
try would split the resonances into two groups, result- 
ing from even and odd /i. The "even" resonances have 
a finite /i = component and disappear upon taking 
the limit R/L — > 0. The "odd" resonances survive in 
this limit, with a finite resonance height and a resonance 
width F ex [R/LY (data not shown). 




FIG. 4: (Color online) Two-terminal conductance of a 
graphene sheet containing a stadium-shaped quantum dot 
without (top) and with (bottom) a 7r-flux. Parameters for 
the calculation are R/L = 0.2, W/L = 12, a/R = ^/3/2, 
d = 2a/3. Without flux tube, the calculation for the conduc- 
tance was done with the method of Ref. 1171 



IV. CONCLUSION 

In this article we investigated the observation of Refs. 
14 16 17, that the two-terminal conductance of a generic 
gate-defined graphene quantum dot shows resonances in 
the limit of a weak coupling to the leads, in spite of the 
naive expectation that electrons can not be confined in 
such a quantum dot because of Klein tunneling. We at- 
tribute this observation to the Berry phase in graphene, 
which quantizes angular momenta to half-integer values. 
With half-integer angular momenta, strict perpendicular 
incidence — the condition for Klein tunneling with unit 
probability — does not occur. As a consequence, con- 
ductance resonances exist in both integrable and chaotic 
geometries. The only difference between the two cases is 
a quantitative one: it concerns the scaling of the reso- 
nance widths with the coupling to the leads.'^ 

The Berry phase can be cancelled against an 
Aharonov-Bohm phase, when a flux tube containing half 
a flux quantum is introduced to the system. With a 
magnetic flux tube, we showed that the relevant angu- 
lar momentum, the kinematical angular momentum, is 
quantized to integer values. In this case a state with 
zero angular momentum is possible. Such a state can 
not form a bound state or give rise to a conductance res- 
onance. We showed this by an explicit calculation for the 
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FIG. 5: (Color online) Behavior of the first three quasi reso- 
nances of the stadium-shaped quantum dot without (top) and 
with (bottom) tt-Aux upon changing the coupling to the leads 
R/L. The other parameters are the same as in Fig. HI 



disc-shaped quantum dot in Sec. [ITj and using numerical 
calculations for disc-shaped and stadium-shaped quan- 
tum dots in Sec. |III[ Once the Aharonov-Bohm phase 
from the 7r-flux tube cancels the Berry phase, the re- 
sults of the full quantum theory are consistent with the 
simple classical expectations. With a tt flux, there is 
a stark qualitative difference between conductance res- 
onance for integrable and non-integrable quantum dots: 
Whereas sharp conductance resonances for the case of 
an integrable quantum dot continue to exist, in the limit 
of weak lead-dot coupling the conductance becomes fea- 
tureless for a generic non-integrable quantum dot. 

In closing, we make two remarks concerning the possi- 
ble realization of the scenario we investigated here. First, 
a flux tube for graphene not necessarily has to be created 
by a real magnetic field, but it can also b e engineered 
via strain as a pseudo-magnetic field.l^ZEll The pseudo- 
magnetic field would have opposite signs for the two val- 
leys, which is of no consequence for our conclusions, be- 
cause the two valleys are decoupled for the smooth con- 
fining potentials we consider here. Second, the appli- 
cation of a well-defined Aharonov- Bohm phase is more 
controlled in ring-shaped structures .^^^ In this case, we 
expect that the our main finding, the strong qualitative 
difference between integrable and non-integrable geome- 
tries in the case of a tt flux, persists. 
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Appendix A: Numerical approach 

The numerical approach follows Ref. [26l The potential 
V{r) is replaced by a potential J^n ^n{y)5{^ ~ ^n) that 
is nonzero at N discrete values —L/2 = xq < xi < X2 < 
. . . < xjv-i < xn = L/2 of the x coordinate, with 

K(y)= / dxV{x,y), n=l,2,...,N-l. 

Between the discrete points the wavefunction is solved 
from the free Dirac equation. This solution takes its sim- 
plest form after Fourier transform with respect to the 
transverse coordinate y, because the free Dirac equa- 
tion does not couple different transverse modes. These 
solutions are then matched by applying the appropri- 
ate boundary conditions at the discrete points x = Xj, 
j = 1,2, . . . ,N. A numerically stable method to imple- 
ment this program is to express both the solution of the 
free Dirac equation and the matching conditions at the 



discrete points x 



in terms of scattering matrices. 



The scattering matrix of the entire sample is then ob- 
tained from convolution of the scattering matrices of the 
27V— 1 individual components. The result of the calcula- 
tion is the transmission matrix t, which is related to the 
two-terminal conductance via the Landauer formula, 



G = 



4e" 



trtff. 



(Al) 



where the trace is taken of the transverse Fourier modes. 
The number of "slices" N and the number of transverse 
modes M must be chosen large enough, that the conduc- 
tance G no longer depends on N and M. 

The vector potential ( 18 1 corresponds to the boundary 
condition 

limip{x, y) = - lim V'(a;, y) for -W/2 <y <0, (A2) 
2;tO xia 

whereas ip{x,y) is continuous at a; = for < y < 
W/2. In the approach described above, this boundary 
condition is expressed in terms of a scattering matrix 
relating incoming and outgoing waves at a; f and x ], 0. 
This scattering matrix has no reflective part, whereas the 
transmission matrix is 



f^ran. — 







-4i/[(fc„-fc„)W^] 



if m - 
if m - 



n even, 
n odd, 



(A3) 



where the integers m and n label the transverse modes 
and kn = 2Tm/W. This transmission matrix has the 
special properties that t = t^ and i^ = 1. 

In order to ensure nmnerical stability, unitarity must 
be preserved while restricting to a finite number of trans- 
verse modes M. For the transmission matrix ( |A3[ ) this 
can be achieved using the following trick: One first 
builds the herniitian matrix h = i(e*'^ — t)/{t + e^"^) — 
cot 4i~t/ sin (j) out of the transmission matrix, where (j) is 



a phase that can be chosen arbitrarily, then truncates 
/i, which can be done straightforwardly without com- 
promising hermiticity, and then uses the inverse relation 
t — e*'^(l -I- ih)/{l — ih) to obtain a finite-dimensional 
transmission matrix. In our numerical calculation we set 
(j) ~ tt/2. We verified that the elements of the resulting 
finite-dimensional transmission matrix approach the ele- 



ments of the exact transmission matrix ( A3 1 in the limit 
that the number of transverse mode Af — >■ oo. 
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